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ABSTRACT 

Deep X-ray surveys by Chandra and XMM-Newton have resolved about 80% of the 2-10 keV cosmic ex- 
tragalactic X-ray background (CXRB) into point sources, the majority of which are obscured AGN. The ob- 
scuration might be connected to processes within the host galaxy, possibly the star-formation rate. Here, we 
use the results of CXRB synthesis calculations as input to detailed Cloudy simulations in order to predict the 
evolution of AGN properties at several mid-IR wavelengths. Computations were performed for three different 
evolutions of the AGN type 2/type 1 ratio between z = and 1: where the ratio increased as (1 + as 
(1 + zf'-^ and one with no redshift evolution. Models were calculated with the inner radius of the absorbing 
gas and dust at 1 pc or at 10 pc. Comparing the results of the calculations to combined X-ray and Spitzer data 
of AGN shows that the predicted spectral energy distributions are a good description of average AGNs found in 
the deep surveys. The existing data indicates that the mid-IR emission from an average AGN is best described 
by models where the attenuating material is ^ 10 pc from the central engine. We present the expected Spitzer 
cumulative number count distributions and the evolution of the total AGN (type 1 + type 2) luminosity function 
(LF) between z — and 1 at rest-frame 8/im and 30/im for the three evolutionary scenarios. The mid-IR AGN 
LF will be an excellent tool to measure the evolution of the covering factor of the gas and dust from z to 1. 
Subject headings: galaxies: active — galaxies: evolution — galaxies: formation — galaxies: Seyfert — in- 
frared: galaxies — X-rays: galaxies 



1. INTRODUCTION 

The earliest optical investigations of active galactic nuclei 
(AGN) identified two major classes of sources: type 1 AGN 
exhibited both broad (FWHM ~ 10^ km s^^) permitted lines 
and narrow (FWHM < lO'^ km s^^) forbidden lines in their 
spectra, whil e type 2 AG Njjre sented only the narrow lines 
iKhachikian & W eedman , 19741) . Unification models of AGN 
propose that this difference can be simply explained by an ob- 
scuring medium that surrounds a basic black hole/accretion 
disk system, and that different lines-of-sight into and through 
this obscuration give rise to the variety of observationjil prop- 
erties witnessed in AGN phenomenology (e.g., lAntonuccil 
119931) . The obscuring medium is often visualized as a geo- 
metrically thick torus comprised of dust and molecular gas 
and of order a parsec in size, sit uated betw een the broad- 
line and narrow-line regions (e.g.. lPier &Kro lik 1992). Ad- 
ditional obscuration is likely to occur in larger structu res as - 
sociated with the host galaxy (e.g.. Maiolino & Riekeiil995h . 
The net unobscured view from the AGN must encompass only 
about 20% of the sky to explain the roughly 4: 1 ratio of ob- 
scured type 2 to unobscured type 1 AGN in th e local Universe 
(e.g., MaioHno & Rieke'1995; Ho et al. 1997). This picture is 
pr obably c onsistent with the recent finding from Sloan data 
by iHao et a l. (2005) that the type 2/type 1 ratio is close to 
unity, because about half of the AGN appear to be well hid- 
den in their host galaxies and would not be seen by Sloan 
but can only be d etected with very dee p spect ra such as those 
used bv lMaioUno&Rieke (19951) and iHo et al. (1997). The 
unification model explains a number of observations. For ex- 
ample, the near infrared (IR) emission of many AGN shows 
a break near 1/xm (e.g., Rieke 1978; Neugebaueret all l 19791 
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lBarvainisiri987t ISanders et al]ll989h . signifying dust that is 
at or near its sublimation temperature and hence close to the 
AGN. Recent Chandra XMM-Newton X-ray observations 
have found that narrow Fe Ka lines are common in both type 
1 and type 2 objects, which implies the presence of cold, low- 
velocity ga s intercepting a significant fraction of the nu clear 
emission dYaaoob & Padmanab han''2004VNandra''2006). 

While the existence of the absorbing gas around AGN 
is well established, the geometry, distribution, origin, and 
evolution of the material is largely unknown. Observations 
of the soft X-ray absorption can measure the equivalent 
neutral column density of the absorbing gas along the line 
of sight, N- H- There is a good correspondence (at the ^ 
80% level; iTozzi et a l. 2006) between the sources with A'h 
> 10^^ cm~^ and those optically identified as being of type 
2. This correspondence implies that, on average, the X-ray 
absorbing gas does lie outside the broad-line region. How- 
ever, observ ations of iVn variabili ty on week- to year-long 
timescales (Risaliti et al. 2002, 2005) argue against a simple 
geometry. Further evidence for a complicated geometry is 
found in the mid-IR. The IR spectra of most AGNs (in vf^ 
units) peak between 3 and 10/im, and roll over at ^ 100/im 
( Andreani et al. 2003; Kuraszkiewicz et al. 2003). This be- 
havior can be described by thermal emission by dust ranging 
in temperature from ^ 40 K to the sublimation temperature of 
1000-2000 K. Part of the cooler emission may be heated by 
star formation rather than the AGN ('e.g.. lMaiolino & Riekd 
1995; Schweitzer et al. 2006). Interestingly, the spectral en- 
ergy distributions of type 1 and type 2 AGN are similar in the 
mid-IR (A > S^m) (Soinoglio & Malkan 1989; Faddaetafl 
[^98; Kuraszkiewicz et al. 2003; Alonso-Herrero et al. 2003 
Eutz et al. 2004; Rigbv et al. 2004), in conflict with the pre- 
dictions of many compact torus models (Pier & Krolik 199^ 
fciranato & Danese 1994; Efstathiou & Rowan-Robinso3 
11995; Granato et al. 1997). All of these observations argue 
for a more distributed or complex geometry to the absorbing 
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material around the AGN (iKuraszkiewicz et al.l2003h . 

If the obscuring material in AGN is distributed in a complex 
manner, then large object-to-object variations in properties 
would be expected, complicating the interpretation of both the 
AGN characteristics and the relationship to the host galaxy (if 
any). Progress can be made, however, by determining the ratio 
of type 2 to type 1 objects (hereafter denoted as R) at differing 
AGN luminosities and redshifts. Under the unified models, 
R is related to the covering factor of the absorbing gas and 
dust, so changes in R give clues to the nature of the obscuring 
structures. For instance, the simplest unified model, that of an 
unevolving torus, predicts no evolution of R. However, over 
the last few years, surveys at multiple wavelengths have noted 
that the fraction of type 1 objects increases with AGN lumi- 
nosity (Uedaetal. 2003; Bareeretal. 2005; Simpson" '2005; 

tio et al. 2005). This trend may simply be a selection effect 
reister et al. 2004), because it is easier to find type 1 high 
luminosity AGN (i.e., quasars) with their broad optical lines 
than type 2 objects, which ma ny models predict to be deeply 
buried (e.g., Hopkins et al. 2005). Moreover, at high redshifts 
dilution from galaxy light can overwhelm the narrow and rel- 
atively weak lin es of a type 2 AGN, increasing the difficulty 
of identification JMoran et akfeOOZ) . Alternatively, the lumi- 
nosity dependence of R may be real, and imply that some 
part of the obscuring material is close enough to the central 
engine to be affected by high luminosities, perhaps through 
dust sublimation (Lawrence 1991) or radiative acc eleration of 
magnetocentrifugal winds cKonigl & Kartjell994) . 

There are also hints that R in creases with redshift at 
moderate AGN luminosities. Both 'Bar ger et al.l (12005 ^ and 
|La Franca et al. (2005), in their analyses of deep X-ray sur- 
veys, found indications for an increase in the type 2 fraction 
from z = to z ^ 1. The sources discovered in these ob- 
servations account for a large fraction of th e cosmic X-ray 
background (CXRB) in the 2-10 keV range ( Mushotzk v et all 
UOOO; Brandt & Hasinger 2005; Worsley et al. 2005), whose 
spectrum can be fitted only if i? ~ 3 — 4 tGillk2004<). 

The implication of these observations is that if the ob- 
scuring medium around an AGN evolves with redshift then 
it must be connected in some way to the evolution of the 
host galaxy. The redshift distribution of the CXRB sources 
peaks at z 1 dTozzi et al. 2001; Bargeret al. 2002, 20051) . 
interestingly close to where the cosmic star-formation rate 
(SFR) density also reaches its maximum (Hopkins 2004). 
This coincidence in redshift distribution suggests a possible 
connection between the absorbing gas and dust around an 
AGN and the SFR wit hin the host g alaxy (already seen in the 
Seyfert 1 NGC 3227:'Davie s et al.l2006 ). With this motiva- 
tion,|Ballantvne et al. ( 2006) determined that a R that evolves 
with both luminosity and z (constrained by the type 1 frac- 
tions observed by Bargeret al. 2005) can fit the observed 
CXRB spectrum and X-ray number counts. However, the 
CXRB itself provides little constraint on the evolution as both 
rapid, slow or zero (e.g., Treister & Urrv 2005) evolution can 
account for the data. A measurement of the true extent of 
any R evolution will provide vital information for models of 
galaxy and black hole evolution over this redshift range. It 
is therefore important to search for other observational sig- 
natures that can differentiate among various evolutions of the 
AGN obscuring material. 

The natural next step is to consider the evolution of AGN, 
in particular obscured AGN, in the mid-IR where the ab- 
sorbed radiation will be re-emitted. With the recent launch 
of the Spitzer Space Telescope, deep sensitive surveys of 



AGN in the mid-IR are now being pe rformed, and are 
finding many high-z (z > 0.5) A GN (fLacv et alJ | 2004| 
iRigbv et al."2004l'Franceschini et alJl2005|ISt ern et al."2005t 
Alonso-Herrero et al. 2006; Barmb vet alJl20 06; Brand et al] 
2006). Therefore, this paper presents predictions for the mid- 
IR number counts and luminosity functions (LFs) for three 
different evolutions of the type 2/type 1 AGN fraction be- 
tween z = and 1. All three evolutions are chosen to fit the 
CXRB, and therefore the mid-IR predictions are consistent 
with the current best X-ray constraints on AGN propert ies and 
evolu tion. In contrast to most previous work (e.g. Silva et alJ 
12004ft . the IR predictions are not based on observed spectral 
energy distributions (SEDs) but are computed using the pho- 
toionization code Cloudy. This method has the advantages 
of being able to directly relate the X-ray properties to other 
wavelengths and to allow interesting physical properties of 
the obscuring medium to be varied and possibly constrained. 

We begin in the next section by reviewing the CXRB results 
of Ballantvne et al. (2006), in particular the three evolutions 
of R that will be investigated here. The following section 
(§13 describes the Cloudy models in detail, including a de- 
scription of our assumptions, and how we average individual 
spectra together before estimating observables. Sect. |3 then 
checks how well these model SEDs compare against the ob- 
served AGN X-ray and mid-IR properties. Sect. |5l presents 
the number counts and LFs for the three AGN evolutions and 
compares the predictions to the available data. In Sect. |6lwe 
attempt to take into account the effects of star formation in 
the model spectra and present revised predictions. Finally, we 
discuss our results in Sect.0and summarize in Sect.ls] 

This paper assumes the standard first-year WMAP A- 
dominated cos mology: Hn = 70 k m s^^ Mpc^^, f^A = 0.7, 
and = 0.3 ^ScSiLSEnni. 

2. REVIEW OF CXRB MODELS 

Wlien fitting the CXRB spectrum, a key parameter is the 
ratio of type 2 to type 1 AGN at a given 2-10 keV X-ray lu- 
minosity Lx and z, denoted R{Lx, z). This can be directly 
related to the fraction of type 2 AGN, /2, at any {Lx,z) 
by /2 = ^/(l + R)- As is commonplace, we define type 
2 AGN to be those that have absorbing column densities 
-^H > 10^^ cm~^ in the X-ray band. If the unified model is 
correct then /2 is approximately equal to the covering factor 
of the > 10^^ cm~^ gas around the AGN. However, the exact 
distribution of this material is unknown. That is, the fraction 
of AGN observed to have the various obscuring columns, the 
iVn distribution, is unknown except for local bright Seyfert 2s 

(e.g., R isaliti et al 1999ft . Under the unified model, all these 

columns of gas exist close to the black hole, with the observed 
iVn distribution giving the relative covering fraction of each. 
This distribution may also be a function of Lx and z. There- 
fore, constructing a CXRB synthesis model reduces to pro- 
ducing a model for how the X-ray absorption may vary with 
luminosity and redshift. 

Two different assumptions on the A'h distribution were 
used by Ballantvne et al. ( 2006) and we adopt the same distri- 
butions here. Ten values of A^h were considered: log iVn = 
20, 20.5, . . . , 24.0, 24.5^ In the 'simple distribution', any 
type 1 AGN has an equal probability p of being absorbed by a 

^ There do exist examples of completely Compto n thick AGN wit h esti- 
mated columns of logiYn > 25 (e.g., NGC 1068: IMatt etailfTggTl) . The 
X-ray emission from these sources is so reduced that they are not an impor- 
tant contributor to the CXRB, but they could potentially add to the mid-IR 
AGN number counts and luminosity functions. However, test calculations 
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column of log A'h — 20 or log A^h = 21.5. Likewise, a type 
2 AGN had an equal chance of being absorbed by a column 

of logiVn = 22 or logTYn = 24.5: 



20.0,. 
22.0,. 



,21.5 
,24.5 



V 
V 



(1 - /2)/4.0 
/2/6.O 



(1) 



The second assumed A'h distribution was the observed 
iRisaliti et" ( 1999) distribution for type 2 AGN, and the sim- 
ple one for type 1 . R isalitj et al,, (.1999) found that in a care- 
fully selected sample of local Seyfert 2 galaxies, 75% had 
log A'h > 23withhalfbeingCompton-thick(logAfH > 24): 



A'h = 



20.0, ...,21.5 
22.0, ...,23.5 
24.0,24.5 



p = (1 - /2)/4.0 
V = /2/8.O 



(2) 



'Ballantvne et al. (2006) considered many different parame- 
terizations for possible evolutions of R(Lx , z), and presented 
results for two cases that could fit the observed shape of the 
CXRB and X-ray number counts (the results were indepen- 
dent of the A'h distribution assumed). The first evolution had 
an initial z = 0, Lx = lO"'^ ''' erg s^^ type 2/type 1 ratio 
of Rn — 4, comparable to previous optical measurements 
dMaiolino & Riekall995l) . and required only gradual redshift 
evolution: 



f2 = Kil + zf'{\ogLx) 



\-4.8 



(3) 



where if is a constant defined by Ro. This evolution of the 
absorbing gas around AGNs also provided an acceptable fit 

(reduced = 1.3) to the type 1 fractions measured by 

[Saraer et al] (l2005h at different luminosities and redshifts. 
The second parameterization presented by Ball antvne et al.l 
l|2006i) was inspired by the SDSS measurement of i?o = 1 
llHao et al.l2 005). and required much more rapid redshift evo- 
lution in order to fit the CXRB: 



/2 = K{1 



X) 



(4) 



This model did a poor job fitting the Barger et al. type 1 frac- 
tions (reduced — 2.6). In both these cases, the z evolution 
was halted at z = 1, because there was no constraint on /2 at 
higher redshifts. Also, the increase in /2 with z may slow sig- 
nificantly at z 1 if the gas obscuring the AGN is connected 
to the cosmic SFR density. 

The CXRB spectrum and X-ray number-counts can also be 
fit with a model where the type 2/tvpe 1 AGN ratio does not 
evolve with redshift (e.g.. iTreister & Ur rv 2005). This sce- 
nario should be considered as a null-hypothesis in our in- 
vestigation. One of the unevolving parameterizations that 
[Ballantvne et al. (2006) presented that fits the CXRB spec- 
trum is 

/2=ifcos — y (5) 

where Rq = 4. This model does have a modest luminosity 
dependence, but this will not affect the comparison between 
the z-dependent evolutions above. 

Figure [2 shows contours of R{Lx, z) for the three evolu- 
tions presented here. The different models have very different 
predictions for the ratio of obscured to unobscured AGN par- 
ticularly at 0.5 < z < 1. The rapidly evolving model (Eq.0} 

showed that the inclusion of such sources did not change our A'^jj -averaged 
SEDs in the mid-IR from which we make our predictions. Thus, objects with 
columns of log A^h > 25 can be safely excluded from our calculations. See 
Sect. l3.5l for details on the construction of the Afn-averaged SEDs. 



predicts ratios of 4: 1 to 5: 1 for quasars at z ^ 1, as opposed to 
the 2:1 to 1:1 ratios predicted by the other two evolutions. At 
lower luminosities and redshifts, the slower evolution model 
(Eq. |3} predicts a larger value of R than the other possibili- 
ties. This paper will explore how these differences manifest 
themselves in the IR AGN number counts and LFs. 

3. CALCULATION OF AGN SPECTRAL ENERGY 
DISTRIBUTIONS 

3.1. General Methodology 

If the gas and dust causing the X-ray absorption are also re- 
sponsible for the majority of the IR emission from an AGN, 
then it may be possible to discriminate among the three evo- 
lutions by mid-IR measurements. In this case, the IR emis- 
sion for a given X-ray A'h can be predicted by a photoion- 
ization model, and since the evolution of A'h is predicted by 
the CXRB models, the subsequent evolution in different IR 
bands can be predicted. The photoionization models also al- 
low the ability to vary parameters such as the location, density 
distribution and dust content of the absorbing gas, which may 
result in constraints to some of these physical properties. 

We employ the ph otoionization code Cloudy v. 05.07.06 
llFerland et alJ Il998h to compute the IR emission of the 
absorbing gas. Although the last decade has seen sig- 
nificant work on the exact radiative transfer problem of 
IR photons through a dust tor us (,Gran ato & Danese 1994 
Efstathiou & Rowan-Robinson | 1995t iG ranato et al. 199^ 
iPuUemond & van BemmeL .20051 iFritz et al.. .2006i) . these 
techniques do not provide an easy correspondence to the X- 
ray absorbing gas. Furthermore, they ignore gas-grain interac- 
tions such as gas-heating by photoelectric emission and grain 
charging by electron capture ( Weingartner & Draine 2001). 
These effects are included in the latest versions of Cloudy 
which also uses a new, more physical grain model where the 
size distributions of both the silicate and graphite components 
as well as RA^Hs are each resolved into 10 size bins. The dust 
module in Cloudy also takes into account the non-equilibrium 
temperature spiking process that affects the very small grains 
(Purcell 1976; Sellgren 1984; van Hoof etal. 2004). The 
computed SEDs therefore include all the detailed dust physics 
of the latest torus models, but also allow a connection to the 
X-ray properties. 

The disadvantage to this technique is that it re- 
quires a simplification of the actual IR emission region 
which may be quite complex (e.g., Nenkova et al.l l2002t 
Dullemond & van Bemmel 2005). Furthermore, there are 
many instances of local AG N whose X-ray absorbing ga s 
seems to be mostly dust-free (IWeingartner & Murravll20()2l) . 
or where only a frac tion of the total column is dusty 
dBallantvne et al.' 2003). In addition, the radiative transfer 
techniques employed in Cloudy are less sophisticated than 
those used in the multi-dimensional torus models (Ferlan^ 
l2003h . However, our goal is not to reproduce the SEDs of in- 
dividual objects, but to make predictions for the average prop- 
erties of an ensemble of AGN at a given Lx and z. Since 
even jn local AGN, there is a wide range of IR properties 
JWeedman e t al. 2005; Buchanan et al. 2006), it is plausible 
that no systematic errors will result from the simplicity of the 
model calculations. This expectation is tested against real data 
in Sect. El 

3.2. Cloudy Model Setup and Assumptions 

Each Cloudy model has a simple setup: a constant AGN 
spectrum with a 2-10 keV luminosity Lx strikes a cloud with 
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an inner radius r pc from the continuum source. The gas 
has abundances similar to the Orion nebula (using the 'abun- 
dances Orion' command in Cloudy), and has a uniform hy- 
drogen density of lO** cm^'^, which is typical for a molecu- 
lar cloud. The calculation proceeded through the cloud un- 
til the X-ray measured column density reached a value iVn 
(utilizing the 'stop effective column density' command'* ). 
The ionizing spectrum was defined using the 'agn' contin- 
uum in Cloudy with an X-ray photon-index of 1 .9, a big blue 
bump temperature of 1.4 x 10^ K (appropriate for 10^ 
black hole accreting at one-tenth of its Eddington rate), and 
an ttox = —1.4. While this is a comm on value of any over 
a wide range of luminosity and redshift, ISteffen et al.l (12006 *) 
showed that this parameter is anti-correlated with the UV lu- 
minosity of the AGN. There is significant scatter in this rela- 
tion, however, and to date only optically selected, very high 
luminosity, z > 4 quasars have been fo und to typically exhib it 
significantly lower values of ctox (e.g.. lShemmer et al.l2006 ). 
Unfortunately, it is unknown if there is any relation between 
ttox and the infrared properties of AGN. 

Cloudy has two stored grain (graphite and silicate) distri- 
butions. One, denoted 'Orion', is based on the flatter redden- 
ing observed in star-forming regions due to a lack of small 
grains. The other, called 'ISM', reproduces the standard 
Galactic reddening law. Since it has been argued that small 
dust grains close to an AGN will be preferentially destroyed 
(|Maiolino et al. 2001; Weingartner & Murrav 2002), we em- 
ploy the Orion grains in our models; however, two grids with 
the ISM grains were also computed to check the sensitivity 
of the results to this assumption. PAHs were included in all 
calculations. 

Finally, consistent with the unified model, we set the cov- 
ering factor of the illuminated gas and dust equal to /2 when 
> 10^^ cm^^ or 1 — /2 otherwise. This assumption may 
run into trouble at high luminosities if, as has been proposed 
(e.g.,|Sanders et al. 1989), quasars evolve over time from be- 
ing completely buried to unobscured. In this model, the cov- 
ering factor of the obscuring gas and dust need not be related 
to the fraction of type 2 sources at any (Lx,z) pair. How- 
ever, for the time being we will continue to assume this is the 
case and that the unified model is valid over all luminosities. 
Sect. l7.4l tests the validity of this assumption. 

Since the covering factor is variable between z — and 
1 when using the evolutions described by equations |3] and |3 
Cloudy models also had to be calculated as a function of z for 
each Lx and Nn. For completeness, the cosmic microwave 
background with the appropriate temperature and intensity for 
the particular z was added to the illuminating spectrum using 
the 'cmb' command. 

3.3. Choosing the Inner Radius of the Absorbing Gas 

The minimum distance from the radiation source that a dust 
grain can typically survive is the sublimation radius, which is 
a function of both the grain size and composition as well as 
the luminosity of the AGN (Barvainis 1987; Laor & Draine 
For a bolometric AGN luminosity of 10 ° erg s , 

or & Draine ( 1993) show that subliming grains exist from 
~ 0.2 pc (for the largest grains) all the way to 4 pc (for the 
smallest grains) from the central engine. Since the Cloudy 
models resolve the grains in both size and composition, we 

^ We confirmed that the X-ray hardness ratio derived from the Cloudy 
models using this method agrees with ones calculated with a wabs*pow 
model in XSPEC. 



are unable to pick an inner radius for the absorbing material 
that will self-consistently account for subliming grains over 
all grain radii. For example, if we chose an inner radius of 
^ 4 pc so that only the smallest grains are subliming then 
we are neglecting the near-IR emission from the larger grains 
that can exist much closer to the AGN. In contrast, if we pick 
a radius of ^ 0.2 pc, the sublimation radius for the largest 
grains, then the smaller grains would be above their sublima- 
tion temperatures and radiating at inappropriately small wave- 
lengths. Self-consistently accounting for grain destruction is 
a complicated procedure especially if it occurs only over a 
fraction of the cloud depth. Therefore, as a compromise be- 
tween the two extremes, we have selected the inner radius of 
the absorbing cloud to be r = 1 pc. For the Orion (ISM) grain 
distributions, sublimation temperatures for the smallest grains 
begin to be reached at \ogLx = 43.5 (43.25) and no attempt 
is made to correct the results for grains radiating above their 
sublimation temperature. This will have the largest effect on 
SEDs computed at the highest luminosities (log Lx }t 46) at 
which a significant portion of the hot dust emission is shifted 
to shorter wavelengths than observed. These high luminosity 
objects are rare, and as our predictions depend on the AGN lu- 
minosity function, this issue will have only a small effect on 
our predictions from z = 0-1. Moreover, Sect. shows that 
when averaged over the A^h distribution, the computed SEDs 
have very similar IR properties to the ensemble of observed 
AGN. 

Choosing a fixed inner radius for the attenuating cloud is 
a different but equally valid technique as that employed by 
recent models of dust tori, which select the inner radius for 
a sing le grain size and change it with AGN luminosity (e.g., 
iFritz et al. 2006). This latter strategy may be more appropri- 
ate for fitting individual SEDs, but we are most interested in 
determining the properties of average AGNs that are being 
selected in the deep Chandra and Spitzer surveys. Thus, it is 
more appropriate to think of the inner radius of r = 1 pc as 
an average distance or as a radius where there is a significant 
density enhancement in an average AGN. 

In this context, we also calculated models with the inner ra- 
dius at r ~ 10 pc. As mentioned in §[0 there is significant 
observational evidence that dust emission over a range of radii 
is required to fit the observed SEDs of many AGN. In addi- 
tion, if the obscuring material is connected to a galactic-scale 
phenomenon such as a starburst ring, then a significant density 
enhancement at distances of ^ 10 pc may be common. It is 
therefore interesting to check if the ensemble of AGNs shows 
evidence for enhanced density at these distances. Clearly, dust 
emission from ^ 10 pc may be too cool to account for some 
of the observed near-IR AGN properties (sublimation temper- 
atures were reached at \ogLx = 45.5 for the Orion grains), 
but it may be important in the mid-IR. 

3.4. The Model Grids 

Our interest is in predicting the mid-IR properties of AGN 
as a function of the 2-10 keV X-ray luminosity Lx and z 
using the fits to the CXRB as constraints. Thus, we iterate 
from \ogLx = 41.5 to 48 (in steps of 0.25), and z — 
to 5 (in steps of 0.05) as we did in constructing the CXRB 
synthesis model (Ballantvne et al. 2006). For each {\ogLx, 
z) pair we compute using the evolution equations |3] S or |5] 
the value of /2, and therefore the covering factor of absorb- 
ing gas around the black hole (if the unified model is cor- 
rect). A Cloudy model is then computed for each log A'h = 
20.0, 20.5, . . . , 24.0, 24.5 cm"^. In practice, because the red- 
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shift evolution is only constrained from z = to 1, we assume 
that there is no redshift evolution in /2 when z > 1 (see also 
§|2}. Therefore, the photoionization models were calculated 
only up to z = 1, and the z = 1 models are then used at 
z > 1. This resulted in 5670 individual Cloudy models for 
the two redshift-dependent evolutions (eqs.|3land|4} and only 
270 models for the unevolving torus model (eq.|5}- 

Grids were computed for all three evolutions with the inner 
radius of the absorbing cloud at r = 1 pc and at r = 10 pc. 
For ease of reference, we will now refer to results from indi- 
vidual grids using a compact label based on the initial Rq and 
r of the grid (e.g., 4:l/lpc, l:l/10pc, or nozevol/lpc). 

3.5. Constructing the Final Spectral Energy Distributions 

At this point, within every grid, there are 10 Cloudy models 
at each (Lx,z) pair for the different A'h values. The pho- 
toionization models provide several different spectra for each 
computation. The ones of most interest for these purposes 
include the sum of the transmitted continuum and outward 
directed diffuse emission from the illuminated material (cor- 
responding to observing the AGN through the cloud), and the 
reflected emission from the inner surface of the cloud. How- 
ever, under the unified model, since all columns exist around 
the AGN, we would observe reflected emission from all A'h 
including the Compton-thick columns. Evidence for this can 
be found in the X-ray spectra of many Seyfert Is where nar- 
row Fe Ka lines are observ ed from optically-thick rnateria l 
out of the line-of-sight (e.g..' Yaaoob & Padmanabhanl2004l) . 
Thus, to construct the SED for each A'h, a weighted average 
of the reflected emission from all 10 A'h models is added to 
the transmitted-Houtward diffuse spectrum for that particular 
A'h- The weights in the average followed the chosen A'h dis- 
tribution, so that for the 'simple A'h distribution' the reflected 
spectrum for all the A'h > 10^^ cm~^ models had a weight of 
/2/6 (eqn.[3. In this way, we calculate a spectrum for each 
A'h that corresponds to observing the AGN through that col- 
umn but also takes into account the reflected emission from 
the other columns that lie out of the line of sight. We refer to 
these SEDs as the 'unified SEDs'. 

The last step before calculating specific predictions for the 
mid-IR wavebands is to average these unified SEDs together 
to find an average AGN SEE) at this luminosity and red- 
shift. This is the same procedure as is used when constructing 
CXRB synthesis models (e.g., Gil li et al.l fT999. 2001). Two 
different A'h distributions (described in Sect. |2j are used to 
generate this final SED. Figure |2] shows examples from the 
4:l/lpc and 4:l/10pc grids when \ogLx = 43 and z = 0.45 
(the spectra are plotted in the rest frame). At this luminos- 
ity and redshift, eqn.|3]gives a type 2 AGN fraction of 0.754 
(R ~ 3), so the A'n-averaged spectrum will be dominated 
by the 'unified SEDs' with A'h> 10^^ cm^^ -j-jje Risahti et 
al. distribution, with half of the type 2 AGN being Compton- 
thick, only makes a significant difference to the average SED 
at wavelengths > 30 jjm. The X-ray part of the averaged SED 
is only moderately absorbed as is required to fit the CXRB 
(see below). 

An interesting property of the A'n-averaged SED shown in 
the top panel of Fig.|2lis the strength of the silicate emission 
at ~ 10/im even though at large values of A'h silicate ab- 
sorption is commonly observed. This emission, as well as the 
majority of the other emission lines, is a result of including 
the reflected emission from the hot inner walls of the absorb- 
ing cloud. The fact that silicate emission is predicted even 
for spectra which are dominated by type 2 AGN is interesting 



given that the majority of AG N with typical S eyfert 2 levels of 
A'h show silicate absorption iShi et al.ll2006l). although some 
exceptional cases do show emission (■Sturm et alJl2006l) . As 
our main focus is on the overall continuum properties, the be- 
havior of this feature, in addition to the emission lines, will 
not significantly influence the predictions. We do note, how- 
ever, that the A'n-averaged SED from the r — 10 pc model 
(lower panel in Fig. |3 exhibits significantly weaker silicate 
emission than in the case where the inner radius of the atten- 
uating material lies at 1 pc. 

4. PROPERTIES OF THE MODEL SEDS 

Given the simplifying assumptions made in these calcula- 
tions, it is important to verify that the A^n-averaged SEDs 
are fair representations of reality, particularly in the mid-IR 
where the predictions of the number-counts and LFs will be 
made. Our goal is not to fit the SEDs of individual objects, 
but to compare various properties of the ensemble of mod- 
els to the large samples of AGN being studied in the deep, 
multi-wavelength surveys. This is valid because we will only 
be making predictions (cumulative number count distributions 
and luminosity functions) that will be compared against data 
that are averaged over many individual objects. 

The photoionization models greatly ease the ability to make 
the comparisons to data, as we are able to calculate the ob- 
served X-ray flux (in addition to the IR fluxes) in multiple 
bands. In this way we can make true multiwavelength com- 
parisons to data. 

The 'simple A'h distribution' (eqn. ^ was used for the 
models presented throughout this section. Employing the 
Risaliti et al. (1999) distribution (eqn. |2j only makes a very 
minor difference to observables like colors, so these results 
are not presented here. 

In this section and throughout the paper, model IR fluxes 
were computed using the predicted flux at the specific wave- 
length of interest (e.g., 3.6/im, 8/xm, etc.). Spitzer of course 
measures an average flux over a passband. Comparing fluxes 
computed at one wavelength with those averaged over the 
Spitzer passband yielded differences of a few percent. Given 
the level of the systematic assumptions made in the calcula- 
tions, this difference in predicted fluxes is negligible. For this 
reason we also do not consider the wavelength dependent ef- 
ficiency of the Spitzer passbands. 

4. 1 . Mid-IR to X-ray Flux and Luminosity Ratio 

If the attenuating gas and dust around an AGN was strictly 
responsible for the thermal IR emission, then a correlation 
might be expected between the relative strength of the IR flux 
and the absorbing column density. That is, the most obscured 
sources should produce the largest mid-IR to X-ray flux ratio. 
However, mid-IR observations of X-ray selected AGN have 
shown no correlation at all between the X-ray hardness ratio 
(a proxy for the absorbing column A'h since more heavily ab- 
sorbed objects will appear hard er) and the strength of the IR 
flux r elative to the hard X-ray (iLutz et alJl2004t iRigbv et alJ 

In Figure |3] we test the existence of this expected correla- 
tion with the 'unified SEDs' computed from the Cloudy mod- 
els. Panel (a) plots the /24/im//2-8 kcV ratio as a function of 
the /2-8 kcv//o.5-2 kcV ratio for all the 'unified SEDs' with 
z < 3 in the 4: 1/lpc grid (the results are not significantly dif- 
ferent using other evolutions). Panel (b) shows the same data 
except only includes the results for models with log Lx < 43. 
The black region shows where models with log A'h < 23 fall. 
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the red area plots models with 23 < log A^h < 24, and the 
blue region denotes models with logA^n > 24. The green 
squares are combined Chandra and Spitzer data on AGN in 
the Chandra Deep Field South (CDFS; Riabv et a lJ2004l) . 

The first important result to note is that only objects with 
columns logA^n > 24 produce the largest relative 24um 
fluxes . This is in agreement with the results of Rig bv et al.l 
i2004h who employed local observed SED templates to pre- 
dict flie /24/xm/ /2-8 kov ratio. Panel (b) of Fig. |3lshows that 
the lowest /24/im//2-8 koV ratios are found from objects with 
\ogLx > 43 if the inner radius of the absorber is at 1 pc. 
However, if r = 10 pc then the absorber is less ionized for a 
given \ogLx, and larger hardness ratios are possible, which 
enable these models to fill much more of the observed pa- 
rameter space (panels (c) and (d)). Interestingly, panel (d) 
shows that when the inner radius of the attenuating material 
is at 10 pc, the lowest /24/xm//2-8 koV ratios are found from 
objects with log Lx < 43, the opposite of what was found 
when r = 1 pc. This can be understood by using Fig. Has 
a guide. The ?- = 1 pc models produce significant hot dust 
emission that peaks at wavelengths < 24/im, while the 10 pc 
models have much cooler dust emission peaking at 24fim or 
even longer (the 2-8 keV flux is basically identical between 
the two cases). Therefore, as the spectra are redshifted, the 
hot dust emission in the r = 1 pc models moves through 
the 24/im band and keeps the mid-IR/X-ray ratio high. In 
contrast, the flux at 24/im would be lower at the same Lx 
and z for the r = 10 pc models since the peak mid-IR emis- 
sion would be redshifted out of this band. The CDFS survey, 
being a pencil-beam deep field, is dominated by AGN with 
\ogLx < 44 (Zhen2 et al. 2004). This would seem to indi- 
cate that the observed-frame 24/im emission may originate, 
on average, from material a few pc removed from the AGN. 

It is also interesting to note the few objects with enhanced 
/24/im/./2-8 kcV ratios that are found at very soft hardness 
ratios. None of the computed 'unified SEDs' were able to ac- 
count for these sources. They may have increased 24/im flux 
due to star-formation heating. However, given the significant 
number of assumptions made in these calculations on ioniz- 
ing spectrum, chemical composition, geometry, etc., it is more 
important that the models account for the broad properties of 
the ensemble. 

Figure 0] plots the same quantities as Fig.|3la) except the 
measurements are made from A'n-weighted spectra. As al- 
luded to in the previous section the X-ray hardness ratios of 
these SEDs do not span a wide range. This is necessary to 
fit the CXRB which is well described by a F = 1-4 power- 
law between 1 a nd 20 ke V (iKushino et alJl2002t iLumb etaP 
120021 iDe Luca & Molend'il l2004 . InFig.l4lwe denote the 
/2-8 kov//o.5-2 kcV ratio for such a power-law with the ver- 
tical dotted line. The red region denotes the /24/im/ /2-8 kcV 
ratio measured from models with z < 1 and log Lx < 44. 
The recent deep X-ray surveys all show that these AGN domi- 
nate die production of the 2-10 keV CXRB JUeda etalJ2003[ 
iBarger etal.l l2005t iLa Franca et"an 120051) . Indeed, the red 
area is bisected by the hardness-ratio appropriate for produc- 
ing the correct shape of the CXRB. 

Fig Elshows that the model SEDs, which were constructed 
based on the assumptions of the unified model, can fully ac- 
count for the distribution of data points from the CDFS. In 
fact, no correlation between /24/im //2- 8 koV and the hard- 
ness ratio is predicted by these models. This is, in part, due to 
the poor translation from hardness ratio to absorbing A^h col- 



umn for surveys that span a large range in redshift. For a given 
A'h, the observed hardness ratio will decrease with incre asing 
redshift. T o correct for this problem. 'Rigbv et"al] J2006 ) and 
iLutz et ai] f2004) plotted the ratio of the absorption-corrected 
2-10 keV luminosity to the 6/im luminosity versus the in- 
ferred A^H for AGN with known or estimated redshifts. These 
authors again found no correlation of the X-ray-to-mid-IR lu- 
minosity with absorbing column at z or at z = 1. Fig- 
ure |5] shows the predicted Lx /i^L^(5.7 fim) versus logA^n 
plot for the unified SEDs resulting from the 4:l/lpc and 
4:l/10pc grids. As the ratios are constructed using rest-frame 
quantities, only one redshift is plotted (z — 0.7). The pre- 
dicted ratios are color-coded based on X-ray luminosity with 
black triangles corresponding to 46 < log Lx < 48, red trian- 
gles denoting 44 < log Lx < 46, and blue triangles showing 
the results for 41.5 < \ogLx < 44. The green squares and 
cyan stars plot th e ratios observed bv lRigbv et alJ J2006t and 
iLutz et alJ (120041) . respectively. 

Both the 1 pc and 10 pc models cover roughly the same 
region in the Lx to vL,A5.1 ^J,m) ratio, but can only explain 
about the upper half of the observed data points. Evidently, 
these unified SEDs do not always produce enough 5.7-6/im 
emission to account for the observatio ns. Star-formation is a 
possible culprit, but tLutz et al.l (120041) correct for this by de- 
composing their spectra using the PAH emission as a guide to 
the star-formation strength. This effect may also be a result 
of omitting the host galaxy from the calculated SEDs, which 
will become more import ant a t low AGN lu minosities (we re- 
turn to this point in Sect. 14. 3> . iLutz et aP (|2Q04) do suggest 
this may be a possible source of error in their measured 6/im 
fluxes, but it is unlikely to account for the roughly factor of 
10 disagreement for some of the points. Perhaps the disagree- 
ment is a result of missing hot dust emission from within 1 pc. 
Replotting the figure with results using the ISM grain distri- 
bution that includes smaller (and hence hotter) grains results 
in only a very small decrease in the Lx to vL^iS.! ^im) ra- 
tio. This result indicates that more hot dust emission may be 
helpful but is again unlikely to solely account for the large 
discrepancies implied by some of the data points. It seems 
more likely that the disagreement is due to a combination of 
all three effects, as well as possible measurement errors and 
the failure of one of our fundamental assumptions in the ge- 
ometry of the absorber A full investigation of this problem is 
beyond the scope of this paper and is left for future work. 

Finally, it is interesting that in the 4:l/lpc models, the 
largest Lx to vLy{5.1 jim) ratios are found at the highest Lx, 
but, for the 4:l/10pc models, they are found at the lowest 
Lx- Indeed, when r — 10 pc, roughly the same spread of 
Lx/i^iy(5.7/im) is found for log > 46 and log < 44. 
When the absorbing material is closer to the ionizing source, 
it is easier to produce hot dust emission at low X-ray lumi- 
nosities, so as Lx increases, there is a slower corresponding 
increase in vL^{5.1 pm). At larger distances, it is more diffi- 
cult to produce hot dust and, for low Lx, increasing Lx also 
results in little response at vL^{5.1 jjm). The sources plotted 
by Rigbv et al. (2006) have \ogLx ^ 44.7, and are thus con- 
sistent with either value of r. 

4.2. Infrared Colors 

Another important check on our calculation strategy is the 
distribution of the observed-frame IR colors. In Figs.|6land0 
color-color diagrams are constructed from the A^n-weighted 
SEDs and compared with Spitzer data presented by Lacv et aP 
(Hool. In the first case, which examines only IRAC colors. 
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the vast majori ty of the models fa ll within the selection region 
determined bv iLacv etal] (12004 . Fig. |7l compares a MIPS 
color to an IRAC one. In both cases, although there are a 
few data points that are significantly redder or bluer than the 
predicted colors, most of the iVn-averaged SEDs have colors 
that are very similar to the observed data. 

The curves traced out by the models in both the color-color 
plots are the changes caused by the IR slope (which depends 
on Lx) being redshifted through the observed bands. For ex- 
ample, the 'tail' at constant /s.T^m/Za.e/im in the l:l/10pc 
panel of Fig.0is caused by AGN with \ogLx < 43 over all 
z. At low Lx and an inner radius of 10 pc, there is little hot 
dust emission, so the SED does not begin to rise until ~ 8/im 
(Fig. The 5.7/im-to-3.6/im slope is therefore nearly un- 
changed over the z range considered here, but /24/im/ /s.T^m 
falls with redshift as less emission gets redshifted into the 
24/im band. 

Another useful way of comparing predicted to observed 
mid-IR colors is t he histogram of v f^{2AiJLm) / v f^{?>^m). 
iBrand et'aP (2006) found that the vast majority of X-ray 
sources in the XBootes survey had vf^{2A^m) / vf^{^^m) « 
0. In Figure m histo grams of this ratio are presented 
from the iVn-averaged SEDs for both the nozevol/lpc and 
nozevol/lOpc grids (black lines), as well as the l:l/lpc and 
l:l/10pc grids (red lines). To more easily compare with the 
[Brand et al. (2()06J result, the histograms only included SEDs 
with /0.5-7 kev > 7.8 X 10~^^ erg cm"^ s^\ the flux limit 
of the XBootes survey, and /^(24/im) > 0.3 mJy. In contrast 
to a strong peak at v Jy{2A^m) / v fu{8^m) w 0, the models 
computed with r = 1 pc tend to have bluer mid-IR colors with 
thepeakat i^/y(24/im)/i^/y(8/^m) « —0.6. The situation im- 
proves slightly when r = 10 pc. In this case, the histogram is 
much flatter with only a slight enhancement of sources with 
negative values of i/f^{24fim)/i/f^{8iim). These results indi- 
cate that models with the inner radius of the absorbing mate- 
rial 1 pc away from the AGN do not produce enough 24/im 
emission to reproduce the mid-IR colors of an average AGN. 
This flux could be added by heating by intense star-formation, 
or, as indicated here, by emission from material at larger dis- 
tances. The 1:1 evolution models produce fewer bluer SEDs 
than the unevolving AGN model. This is likely a result of the 
greater fraction of obscured large-Lx AGN at z ~ 1. 

4.3. Absolute Fluxes 

The final test to be performed here is comparing the ab- 
solute observed-frame fluxes predicted by the A'n-weighted 
SEDs with the available data. Figure |9] plots the predicted 
8/im and 3.6/im fluxes against the predicted 0.5-10 keV 
fluxes. The black crosses are from the 4:l/lpc models, while 
the blue triangles plot the predictions from the 4: 1/lOpc mod- 
els. Only models with z < 3 were plotted on the Figure. To 
compare with these predictions, the IR and X-ray data from 
sources in the Extended Groth Strip (EGS) were taken from 
[Barmbv et al. (2006) and are plotted as red squares. 

Fig-Elreveals an important limitation in our A'n-weighted 
SEDs. As expected, there is a strong correlation between the 
IR and X-ray fluxes in our model, and this does a good job 
describing the distribution of data from the EGS at high X- 
ray fluxes. At low X-ray fluxes, the predicted IR flux severely 
underestimates the average measured fluxes, in particular at 
3.6/im. This is because we do not include any emission from 
the host galaxy in our model SEDs. This will become im- 
portant at the lowest AGN luminosities, especially at 3.6/im 
where there will still be an important contribution from the old 



stellar population of the galaxy. This problem does not seem 
to impact the 8/im fluxes as strongly, but we will not be able 
to make accurate predictions at 3.6/xm or 4.5/im. As seen in 
the following sections, this will not be a significant hindrance 
in producing strong predictions at longer wavelengths. 

As evidence for this, Figure fTolp lots the predicted 2-10 keV 
fluxes from the z < 3 models in the 4:l/lpc and 4:l/10pc 
grids against the 24/im flux. The symbols are the same as 
the previous plot exc ept the red squares plot da ta from the 
CDFS presented by Alon so-Herrero et a l. (2006). The locus 
of model predictions bisects the data points indicating that 
the calculated 24/im flux gives a good description of an av- 
erage AGN seen in the deep surveys. Also plotted in the fig- 
ure are regions where bright local AGN ('Piccionotti AGN') 
and local starburst galaxies are expected to be found (see 
Alonso-Herrero et al. 2004, 2006). The TVn-weighted SEDs 
predict fluxes that place them within or on the lower boundary 
of the local AGN properties, again indicating that the models 
are good descriptions of AGN over a range of redshifts and 
luminosities. 

In summary, the computed SEDs seem to adequately de- 
scribe average AGN behavior from the X-ray to the mid-IR, 
a property that has previously only been possible using SED 
templates. These results also indicate that the assumption of 
a constant aox is not resulting in a serious bias or error Our 
computational approach has the advantage of being able to 
vary the physical properties of the absorbing gas to see which 
arrangement best describes the average AGN. The successful 
tests presented in this section validate this strategy and give 
confidence that the calculated number counts and LFs will 
have significant predictive power. 

5. PREDICTIONS FOR BARE AGNS 

This section presents the predicted cumulative number 
count distributions and luminosity functions resulting from 
the final A^H-averaged SEDs. Recall that our goal is to use 
the results of the deep hard X-ray surveys and the CXRB to 
predict the IR properties. Therefore, a key ingredient in com- 
puting these integrated quantities is the hard X-ray luminosity 
function (HXLF) and its evolution with z and Lx- As before 
(Ballantvne et al. 2006), we employ the luminosity-dependent 
density evolution HXLF of Ueda et al. (2003). This type of 
evolution has also been inferred bv Hasinger et al.. (.2005.) and 
iLa Franca et al.l (l2005h using other X-ray surveys. 

5.1. Number Counts 

The expression used to calculate the cumulative number 
counts at differe nt IR wavele ngths is the same as equation 3 
from Ballantvne et al.l (l2006h . However, an important differ- 
ence in the calculation is that the integral over X-ray lumi- 
nosity is begun at the minimum Lx that results in an IR flux 
S. In this way, for each IR wavelength of interest, number 
count distributions are computed for each of the three R evo- 
lutions considered, and for models with inner radii at 1 and 
10 pc. Only the results with the simple A'h distribution are 
shown below , as the re were only minor differences when the 
Risaliti et al. '(^^9?) distribution was assumed. 

Figureirriplots the predicted cumulative number counts dis- 
tributions at 5.7, 8, 24 and 70/im. The black lines denote the 
three evolutions computed with the inner radius of the absorb- 
ing gas and dust at r = 1 pc from the ionizing source, while 
the red lines show the results when r = 10 pc. Within each 
color, the three line styles differentiate among the three evo- 
lutions of R. The data points in the 5.7, 8 and 24/xm panels 
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are the measured number counts of X-ray selected AGN from 
the GOODS survey ( Tre ister et alJ20 06). 

In the 5.7 and 8/im panels there are two dashed 
green lines that plot the computed number counts for the 
nozevol/ISM/lpc and nozevol/ISM/lOpc models. The smaller 
grains included in the ISM model enhance the short wave- 
length IR emission, and therefore slightly increase the num- 
ber counts at a given flux for these two Spitzer bands. This 
difference is small enough that it will not affect the arguments 
presented below. 

It is clear from this figure that, apart from a small amount at 
24/im, there is little to no difference in the predicted number 
count distribution for the three evolutions of the type 2/type 1 
ratio. This is a result of integrating over Lx and z which will 
tend to wash out any changes among the evolutionary mod- 
els. Differences in the number count distributions do arise at 
lOfim, but they are small in magnitude and it may be diffi- 
cult for future surveys to reduce the errorbars to small enough 
values to discriminate among the different AGN evolutions. 

The number count distributions do show a significant dif- 
ference between the ones based on calculations with ?- = 1 pc 
(black lines) and the models with r = 10 pc (red lines). Be- 
fore drawing conclusions it is important to consider the conse- 
quences of both our methodology and the observatio nal se lec- 
tion effects on this plot. As mentioned above (Sect. 14. 3> . the 
omission of starlight from the host galaxy results in the mod- 
els underpredicting the near-IR flux at low AGN fluxes. This 
will result in an underprediction of the number counts at small 
IR fluxes, with the difference increasing both in strength and 
in flux as we move to shorter wavelengths. The GOODS sur- 
vey covers a small angle on the sky and the data plotted here 
are X-ray selected, so the survey will be missing both bright, 
rare objects and heavily obscured AGN, both of which are in- 
cluded in our calcula t ions. Indeed, Donlev et al. (2005) and 
lAlonso-Herrero et al .! (2006) both find that a significant frac- 
tion (~ 50%) of AGN have been missed in the X-ray surveys, 
even at exposures > 1 Ms (see also iBallantvne e t al. 2006). 
The addition of rare, bright objects will result in the models 
overpredicting the number counts at the bright end. Compton- 
thick AGN, while very faint in the 2-8 keV X-ray band, can 
produce very typical AGN fluxes in the IR (Fritz et al. 2006). 
The inclusion of these sources in our models will result in an 
overprediction compared to the data points at all fluxes, where 
the amplitude of the error is related to the relative number of 
Compton-thick AGN. 

With the above points in mind. Fig. II II indicates that the 
r = 10 pc models provide the best description of the observed 
number counts at all wavelengths. This is perhaps not surpris- 
ing at observed-frame 24/im as it is consistent with our earlier 
conclusions from Fig. |S] that the r = 1 pc models do not on 
average produce enough 24/im flux. However, this conclusion 
at observed-frame 8/im seems at odds with Fig.|9]which indi- 
cates that the r = 1 pc models are a better description of the 
8/fm flux. With this interpretation, the number counts should 
be observed to be underpredicted by our models at low 8/im 
flux, but this is only the case for the r = 10 pc models. The 
two plots can be brought into agreement only if a non-AGN 
emission component begins to dominate the 8/im flux below 
a certain X-ray flux. Indeed, the 8/im flux of the EGS data 
is nearly constant below /0.5-10 kcV ^ lO^^"' erg c m~^ s~^ 
(this i s perhaps clearer in the original Fig. 7 of Barm bvet alJ 
|2p06). Therefore, the r = 10 pc models are consistent 
with both the EGS 8/im fluxes and the 8/im GOODS num- 
ber counts. 



To summarize this section, we have found that the cumula- 
tive number count distributions are very sensitive to the geom- 
etry of the obscuring material of the AGN, and the GOODS 
data are best described by the models where the attenuating 
material is ^ 10 pc from the AGN. It will be interesting to 
compare these models to future AGN counts at 70/im. 

5.2. Luminosity Functions 

As seen in Fig.Q the three evolutions of R considered here 
predict very different values of the type 2/type 1 ratio with 
redshift. It is therefore expected that the redshift evolution of 
the IR luminosity functions should provide a means to dis- 
criminate among different evolutions of the AGN covering 
factor. 

I iu' Ueda et al. (2003) HXLF tells us how the number den- 
sity of AGN per increment of log Lx changes with z and Lx- 
The following expression is then used to convert this to an IR 
rest-frame LF d^/d{log vL^) at a given wavelength, 

^ djlogLx) 
d(log i^L^) d{\og Lx) d{log vL^) ' 
As our photoionization technique allows a straightforward 
comparison between the hard X-ray and IR luminosities for 
each SED, the conversion factor in Eqn.|6lcan be easily com- 
puted. As an example. Figure El plots the input 2-10 keV 
luminosity versus the IR luminosity for three different wave- 
lengths using the results of the nozevol/lpc and nozevol/lOpc 
grids. The range in IR luminosities is smaller than the range in 
the X-ray band, with the IR range decreasing with increasing 
wavelength. For example, at r = 1 pc, although the X-rays 
span over 6 decades in luminosity, the resulting IR luminosity 
spans only about 5 decades at 15/im and 4 decades at 30/im. 
This wavelength-dependent sensitivity on the input X-ray lu- 
minosity is a result of our specific source geometry and the 
resulting amounts of dust available to emit at the different IR 
energies. Emission at shorter IR wavelengths is dominated 
by hot dust, and as the luminosity of the ionizing source in- 
creases, greater amounts of dust are raised to high tempera- 
ture. The longer wavelengths are dominated by cooler dust, 
but as the luminosity increases our models may truncate the il- 
luminated cloud before large amoun ts of this dust are warmed. 
This behavior is confirmed in Sect. 17.41 where we show that 
a reduction in the density of the IR-emitting region, and an 
accompanying increase in its size, can significantly increase 
the long wavelength output of the model at high Lx- 

The factor d(log Lx ) / ^(log i^Liy) used in eq. |6lis of course 
the slope of the lines plotted in Fig. El At a given z and 
wavelength, the IR LFs are then computed for each evolu- 
tionary model by iterating over Lx, finding the appropriate 
IR luminosity, determining the slope d{\og Lx) / d{\ogiyLi,) 
at that luminosity, and then multiplying the HXLF at (Lx , z) 
by that slope. As with the number count distributions, the re- 
sults described below assume the simple A^h distribution as 
described in Sect.jJl 

The predicted rest- frame 8/im LFs for the three R evolu- 
tions are plotted^ in Figure^l^t z = 0.2,0.6 and 2.0. The 
left- and right-hand panels show the r = 1 and 10 pc re- 
sults, respectively. To emphasize any differences among the 
LFs, the lower window in each panel shows, for each of the 
three redshifts, the ratio between the evolving LFs and non- 
evolving LF. It is clear that the LFs do allow a mechanism to 

^ These LFs and the others at longer wavelengths were computed using the 
results from grids with Orion grains. LFs computed with ISM grains made a 
negligible dilference, even at 8/im, and therefore are not presented. 
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determine among the different evolutionary scenarios for the 
AGN type 2/type 1 ratio. If the inner radius of the attenuat- 
ing gas and dust is on average ~ 1 pc from the AGN, then, at 
z = 0.2, both the 1:1 and 4:1 evolving models predict 8/im 
LPs smaller than the non-evolving case by factors greater than 
2 for objects just past the knee of the LF. At higher z, the large 
fraction of type 2 objects predicted by the 1:1 evolutionary 
model (see Fig.^b)) results in a factor > 3 enhancement of 
the LF over the non-evolving model. If the obscuring material 
has an inner radius of ^ 10 pc then these differences move to 
higher luminosities and, at low z, reduce in amplitude. This is 
because when the gas is further away larger AGN luminosities 
are needed to produce the hot dust emission at 8/im. 

The data points in Fig. 1131 are the m easured type 1 AGN 
LF at z — 2 from Brow n et aP J2006h . These objects were 
observed at 24fim and lie at 1 < z < 5 and therefore corre- 
spond to rest-frame emission at 4-12/im. It is obvious that the 
z = 2 LFs computed with r = 10 pc provide a better match to 
the observed LF than the ones computed with r = 1 pc, con- 
sistent with the conclusion on r reached by considering the 
number counts ("Sect. l5Jl . The close match in amplitude be- 
tween the predicted and measured LFs is probably fortuitous, 
as the computed LF is for both type 1 and type 2 AGNs com- 
puted at one z, while the measurements are for only type Is 
spread over a large range of z. Nevertheless, the agreement in 
the shape of the LF with the r = 10 pc models is significant. 

To explore the behavior of the absorbing material at longer 
wavelengths, we plot in Figure^Jthe predicted LFs for rest- 
frame 30fj,m. The differences among the three evolutions are 
enhanced at this longer wavelength, with both the evolving 
models underpredicting the z = 0.2 non-evolving LF by fac- 
tors > 10. Moreover, for both the r — 1 and 10 pc set of 
models the differences persist over a much larger range of lu- 
minosity than at 8/im. At z = 1, the large number of obscured 
AGN predicted by the 1 : 1 evolution leads to an enhancement 
of the 30/im LF by factors that can approach 50 if r = 1 pc 
or 100 if r = 10 pc. The largest differences among the evo- 
lutions appear above the knee of the LFs and move to higher 
luminosities in the r — 10 pc models. Evidently, the warm 
dust emission is very sensitive to the covering factor of gas 
around the AGN, and therefore provides a clear tracer of any 
evolution with z. 

In order to compare these computed 30/t m LFs to data w e 
made use of the SWIRE survey of Francesc hini et alJ i2005l) . 
These authors presented AGN SEDs from Spitzer IR to op- 
tical and X-ray wavelengths for X-ray detected AGN in the 
ELAIS Nl region. Using their SEDs, we estimated the 60/im 
flux density and computed the rest-frame 30/im LF at z = 1. 
The LF was derived by the 1 /Vinax method. The redshift dis- 
tribution of sources peaks around z ~ 1 and the derived LF 
is not sensitive to the choice of the redshift bin. A redshift 
bin of 0.8-1.2 was selected in order to ensure a reasonable 
number of objects (the LF changes by < 20% for a redshift 
bin 0.9-1.1). The range of 30/im luminosity, i^L^y, of objects 
at 0.8 < z < 1.2 is about one order of magnitude and we 
employed a range of i^L^ = lO^'^-lO'*^'^ erg s^^ to exclude 
an outlier. For simplicity, we assume the same incomplete- 
ness over both the redshift and luminosity bins, although this 
can result in overestimating the LF. However, the LF with- 
out the incompleteness correction is smaller by only 30%, 
within the estimated uncertainty. Corrections were made for 
the 99% completeness of X-ray detection (see 'Manne rs et all 
12003) . the 90% completeness of infrared source detection, 
and the identification of 102/122 X-ray sources in the Spitzer 



IRAC image. The uncertainty of the LF includes the Poisson 
noise statistics on the number of sources and the uncertainty in 
the observed 30/im luminosity estimated by fitting the SEDs. 
While the IR SED fitting is based on the IRAC data, most of 
the objects used for constructing the LF also have 24/im de- 
tections and they are roughly consistent with the SED fitting. 
The /25/xm//i2/xm ratio of local ( z < 0.2) Seyfert galax - 
ies varies only by a factor of ~3.5 dNeugebauer et alJll986^ ■ 
Thus, we argue that this factor is a conservative estimate of the 
uncertainty in the observed 30/im luminosity, because the IR 
SED fitting performed by Franc eschini et al. (2005) is based 
on IR data points at several wavelengths. The final LF is plot- 
ted in Figure ^5 where the uncertainty is dominated by the 
uncertainty in the IR luminosity. 

This datapoint is consistent with the predictions of the 
r = 1 and r = 10 pc models. If the rest-frame 30/im emission 
from an average AGN with this luminosity at z = 1 arises 
from material with an inner radius of 1 pc from the cen- 
tral source, then only the 1:1 evolution predicts enough ab- 
sorbed AGN to account for the observed LF. Converting the 
IR luminosity of this data point to an X-ray luminosity gives 
log Lx ~ 44. At this luminosity and 2 = 1, Fig.^shows that 
evolution|3predicts a type 2/type 1 ratio of ~5:1, as opposed 
to the ^2:1 ratio derived from the other two models. These 
lower ratios are consistent with the data point if the inner ra- 
dius of the absorbing gas and dust is at 10 pc. Unfortunately, 
in this situation the luminosities of the AGN are too low to 
distinguish among the three evolutions. 

Moving to even longer wavelengths further enhances the 
differences among the three different models of R{L, z), al- 
though complications arising from star formation are more 
important at longer wavelengths (see next section). The rest- 
frame 35/im LF, which is accessible to Spitzer at z = 1 us- 
ing the 70/im MIPS band exhibits a similar behavior as seen 
at 30/im with the differences in the LFs increased to factors 
- 100 at both z and 1. 

6. ACCOUNTING FOR STAR FORMATION 

In the previous section we found that large differences 
among the three type 2/type 1 evolutions are found at 30/im 
(for both the 1 pc and 10 pc models). However, these longer 
wavelengths will likely be contaminated with emission from 
star-forming regions. This will be particularly troublesome at 
z ^ 1 if the AGN absorbing gas is connected to the SFR in 
the host galaxy. Therefore, some limits need to be placed on 
the effects of star-formation on the above predictions. 

A simple \yav to proceed is to make use of the work of 
lYun&Cariliilll2()()l . who fit a single temperature dust model 
to the far-IR/sub-mm emission of local starburst galaxies 
(they also included bremsstrahlung and synchrotron emission, 
but emission from these processes is only important at much 
longer wavelengths). The median dust temperature Td and 
emissivity index (3 they find from their sample is Td = 58 K 
and f3 = 1.35. By relating the total 40-500/im emission to 
the SFR by SFR= L40-500/im/(5.8 x IO^Lq) Mq yr-i, 
Yun & Carilli (2002) then write the dust emission as 



= 1.56xl02V3SFR 



(1 



-(i//2000)^ 



(•gO. 000831^ _ 2) 



erg s ^ Hz ^ , 



(7) 

where i/ is in GHz. 

As this spectrum assumes only one relatively cool dust tem- 
perature and ignores more complicated dust physics such as 
PAH emission and non-linear temperature spiking from very- 
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small grains, it will not be an accurate correction to the AGN 
spectra at wavelengths < 20//m. It should be a reasonable 
correction at 30^m. The dust emission spectrum (Eq. 
was added to the final A^h -averaged Cloudy SEDs between 
3 and 1500/im assuming constant SFRs of 0.5, 1, 5, 10 and 
100 Mq yr-\ 

As the number count distributions did not provide a good 
method to discriminate among the three evolutionary models, 
we will omit a discussion on how star-formation impacts the 
number counts. Suffice it to say that the only non-negligible 
effect on the predicted distributions occurs at 70/im. 

Turning to the LPs, Figure [21 plots how one particular 
30/im luminosity function is altered by adding in dust emis- 
sion from the different SFRs. For SFRs < 10 M0 yr~^ sig- 
nificant deviations occur only when \ogvLi, < 44.5. At low 
AGN luminosities the majority of the 30/im emission is from 
the star-formation emission. In contrast, the 30/im LF can 
be dramatically changed at nearly all luminosities if the SFR 
~ 100 Mq yr^^. When the SFR reaches these values, the 
30/im luminosity is only mildly sensitive to the heating from 
the AGN. Although there is significant scatter from object to 
object, such a high SFR seems to be reached in an average 
galaxy only at z > 1.5 iTuneau et alJl200.^ iPapovich et alJ 
I2OO6). However, in spite of the changes to the 30/tm LF 
caused by star-formation. Figure^] shows that the LFs from 
the three AGN type 2/type 1 ratios are still significantly differ- 
ent for SFRs < 10 M0 yr~^, and the discussion from the pre- 
vious section is still valid. This figure does suggest, however, 
that generating a LF from sources with a lower SFR would be 
preferable. 

7. DISCUSSION 

The goal of this paper was to investigate how three possi- 
ble evolutions of the AGN type 2/type 1 ratio (all of which 
fit the CXRB) manifest themselves in the mid-lR. The predic- 
tions were computed by employing the photoionization code 
Cloudy to calculate the SED from an obscured AGN over a 
large range of X-ray luminosity, Lx, and redshift. In this way, 
a direct connection between the observed X-ray and IR fluxes 
could be predicted. The covering factor of the absorbing gas 
and dust changed with z and Lx as predicted by CXRB syn- 
thesis models and the unified model. 

The results presented in the above two sections show that 
this photoionization-based method is very successful. The 
family of computed SEDs indeed capture many of the aver- 
age properties in typical AGN observed in the deep surveys of 
Chandra, Spitzer andXMM-Newton, including infrared colors 
and observed IR and X-ray fluxes. Typically, investigations 
of the IR properties of AG N involve either deta iled dust ra- 
diative transfer models (e.g.-' Granato et alJl997h that provide 
no connection to the X-rays, or are based on phenomenolog- 
ical SEDs and are unable to provide constraints on physical 
parameters (e.g.. §ilva et al. 2004 ). Our approach is not de- 
signed to fit the SEDs of individual objects, but we can use 
it to constrain physically interesting quantities in the ensem- 
ble of objects. In this way, a picture of an average AGN 
can emerge, uncluttered from the significant source-to-source 
variability that affects the study of small samples. 

7. 1 . The Location of the Absorbing Gas 

Two different inner radii were used in computing the 
Cloudy models, 1 pc and 10 pc (see ij l3.3> . as tests for whether 
large samples of AGN show evidence for dust emission at 
these average distances. Although only these two radii were 



considered here, comparing the results from the two sets of 
models with the available data gave rise to the interesting 
result that almost all observational tests preferred the 10 pc 
models over the 1 pc ones (Figs. |3] and |8}. In particular, 
the GOODS number-count distributions were best described 
by the 10 pc models from observed-frame 5.7/im to 24/im 

(Fig. ED- 

It is important to note here that both the r — 1 and 10 pc 
sets of models were assumed to have the same constant den- 
sity of lO"' cm^'^. Thus, all Cloudy runs with logiVH < 23.5 
have pathlengths less than 10 pc (ignoring ionization effects). 
As a result the r = 1 pc models will be dominated by hot 
dust, as the column does not extend to great enough dis- 
tances to produce significant warm emission. A new lower 
density distribution that would have an inner radius of 1 pc 
and be able to produce the correct 24/im flux is certainly 
conceivable. On the other hand, if the obscuring medium 
is clumpy with each clump being roughly the same density 
(Dullemond & van Bemmel 2005), then our setup may not be 
too unrealistic. An investigation using different density distri- 
butions will be the subject of future work. 

Although the 10 pc model seems to be the best description 
of data observed at the mid-lR Spitzer wavelengths, the ob- 
served SEDs of many AGN show that hot dust from inside 
this radius is also required. The 10 pc Cloudy models do pre- 
dict dust temperatures close to and at the sublimation value 
at the inner radius of the cloud, but only for the more lumi- 
nous AGN models (see S I3.3> . Thus, there must be dust emis- 
sion over a range of radii and temperature s, in agreement with 
SED fitting of individual objects (e.g., Kur aszkiewicz et all 
|2003). This conclusion is also consistent with observa- 
tions searching for variability in the IR continua of AGN 
that find little to no variations at mid- to far-lR wavelengths 
on day to month long timescales ( Rieke & Lebofskv 198 ij 
'Edelson & Malkan"198T; "Cle ment et al. 1988). In contrast, 
near-lR (JHK band) variability of local Seyferts is observed 
with measured timelags from the optical/UV placi ng the emit- 
ting region < 1 pc from the continuum source ( Clavel etall 
1989; Suganuma et al. 2006; Minezaki et al. 2006). 

Neugebauer & Matthews ( 1999) report long term monitor- 
ing of nearby, low-luminosity quasars at 10.6/im. At the typ- 
ical redshifts of their sample, this wavelength is very close to 
8/im in the rest-frame. They find slow, multiple-year variabil- 
ity but conclude that in no case are the results in contradic- 
tion to a thermal origin for the emission in the objects with 
luminosities < 10^^ Lq. Specifically, the variability time 
scales are consistent with our conclusion that the 8/im emis- 
sion arises from thermal re-radiation by dust in a component 
with a typical inner radius of 10 pc. They also summarize 
previous measurements of luminous, radio-loud quasars, such 
as 3C 273; in this object in particular, observations at ^ 8/im 
in the rest frame show variations that are too fast to interpret as 
thermal radiation. Since such objects are exceptional among 
the total quasar population, it is unlikely that such behavior 
would strongly affect our results for the ensemble properties 
of these sources. 

In the future, model grids can be extended to even larger 
inner radii to establish a pair of bracketing values for this 
key parameter. It will also be possible to break observed 
data into redshift and luminosity slices to determine if there 
is an evolution of the inne r radius with z and/or Lx (e.g., 
iHatziminaoglou et alJ2006l) . 

7.2. Determining the Evolution of R 
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Sections |5] and |6l showed that the mid-IR LF of AGN will 
be able to discriminate among different models of R, the type 
2/type 1 ratio. Out of the wavelengths considered here, the 
rest frame 30/im LF between z = and 1 gives the largest 
differences among evolutionary scenarios, as well as being 
the least contaminated due to star-formation emission. This 
last point needs to be confirmed with a more sophisticated 
prediction of the IR emission due to star-formation. 

Assuming the 10 pc models are a more accurate guide, 
then the right-hand panel of Fig.^Jindicates that largest dif- 
ferences among evolutions of R occurs at luminosities with 
log vL^ > 45. This is past the knee of the LF, and so will re- 
quire wide surveys to pick up the large numbers of luminous 
AGN. The observations will also require significant depth in 
order to trace the LF out to z 1. 

Another issue in observationally testing the LF predictions 
is accessible wavebands. To really understand the evolution 
of the obscuring region, one would like to trace the LF from 
low redshift (say, z = 0.2) Xo z ^ 1. The LF predictions are 
in the rest-frame, which means that for 30/im observations at 
z ~ 0.6, the required flux is being radiated at 48/im. One 
would then be forced to perform multi-wavelength observa- 
tions and SED fitting for each object in order to predict the 
flux at the required wavelength. This SED fitting could be 
relatively simple as the largest uncertainties for determining 
the LF will be in the space density statistics. 

Given these difficulties, it is natural to ask if the mid-IR LF 
is the easiest method to determine the evolution of R. A per- 
haps more straightforward strategy would be to measure the 
ratio of narrow-line to broad-line AGN in the optical as func- 
tions of z and L^q m], the luminosity in the [O 111] A5007 
line. Alternatively, we could measure the ratio of X-ray AGN 
with absorbing column densities log A^h < 22 to those with 
log A^H > 22 as a function of z and Lx, although an accu- 
rate value of iVn requires an X-ray spectrum with sufficient 
counts. Since it is currently difficult to obtain enough counts 
for spectral fitting from faint sources, hardness ratios (which 
assume a spectral slope) are often used as an alternative. 
Thus, both of these other methods to track the evolution of 
R require serious observational investments, but, more prob- 
lematically, they also are both subject to severe selection ef- 
fects. As is becoming increasingly apparent, optical selection 
of AGN (i.e., selection on the basis of high-ioniza tion lines) 
has significant limitations (e.g.- lNetzer et alJl2006l) and many 
(mostly type 2 AGN) are missed due to obscuration or dilu- 
tion ©oran et al. 2002; Eckart et al. 2006; Rigbv et al. 2004). 
Likewise, X-ray selection is limited because the observable 
passband for the most sensitive telescopes ends at ^ 10 keV, 
and the observations are thus insensitive to objects with high 
column densities, especially those with A'h> 10^'' cm^^, the 
Compton-thick AGN. The high energy detectors on Swift and 
Integral have been able to conduct surveys, but are li mited to 
the ne arby Universe due to low sensitivity ( Markwardt et al.l 
l2005t lBassani et al. 2006). Future observations with the Hard 
X-ray Detector on board Suzaku may help, but it is a non- 
imaging instrument. The cancellation of the hard X-ray imag- 
ing mission NuStar makes finding significant numbers of the 
Compton-thick AGN by their hard X-ray emission practically 
impossible until the launch of Constellation-X. 

The advantage of the mid-IR method described here is that 
it is less affected by the above selection effects. The pre- 
dicted LFs shown in Fig. [H] are for both type 2 and type 1 
AGN, therefore detailed optical emission line spectra are not 
required. Compton-thick objects should also be found in the 



infrared, as the absorbed X-ray radiation will be re-emitted 
by the dust. The difficulty, of course, is determining which 
IR sources are predominantly powered by AGN as opposed 
to star-formation, although methods that examine the spec- 
tral shape of the nea r- and mid-IR have had some success 
dAlonsQ-Herrero et al..2006; Brand et al. 2006), and A GN are 
well known to have very weak R\H emission ( Weedma n et alJ 
120051) . However, X-ray selection is still the best way to iden- 
tify AGN, as Compton-thin objects can be observed to high 
redshifts. Deep infrared observations are then required to find 
the Compton-thick sources. Comparisons between multiband 
observations and simulations similar to those performed here 
will be useful in determining the AGN contribution to these 
potential Compton-thick sources. This can be done by adding 
in contributions from the host galaxy to the photoionization 
models, and computing a large grid of IR properties over a 
range of AGN luminosities and redshifts for comparison to 
both near- and mid-IR data. 

Another interesting effect seen in the IR and not in the other 
wavelengths is the de pend ence on the geometry of the obscur- 
ing medium. Figures [T3l and [T4lhave shown different shapes 
for the LF depending on the inner radius of the absorbing gas 
and dust. For the one data point at 30/im, it was difficult to 
determine the best evolution of R without knowing the ge- 
ometry. Fortunately, there are significant differences in the 
predicted LFs and number counts for the different inner radii. 
Therefore, it is possible to break the degeneracy between r 
and R, as was done at rest-frame 8/im (Fig.ll3>. 

Finally, the mid-IR LFs presented above are relatively in- 
sensitive to the actual A^h distribution assumed. We con- 
sidered both a 'simple' A^h distribution, where all column 
densities had equal weig ht within the type 1 or 2 designa- 
tions (eqn.HJ, and the Risaliti et al.1 0199 9) iVn distribution in 
which half of the type 2 AGN are Compton thick and 75% 
had logiVH > 23 (eqn.|2}. These weights were used in con- 
structing both the 'unified' and A^H-averaged SEDs. Despite 
this, the two distributions made only negligible differences 
to the predicted number counts and LFs. This is because 
the additiona l weigh t afforded to Compton thick AGN in the 
'Risaliti et al." ('1999") A'h distribution has only a small effect 
in the mid-IR (Fig.|2j. Therefore, although mid-IR LFs will 
not be able to quantify the relative number of Compton thick 
AGN, they are still very sensitive to the overall evolution of 
the AGN type 2/type 1 ratio. 

7.3. The Role of Star Formation 

Star formation is a key process in this investigation in two 
distinct ways. First, it acts as a contaminant of the mid-IR 
emission from the obscured AGN, masking the evolutionary 
effects we are trying to investigate. Using a straightforward 
estimate of the effects of star formation on the mid-IR, we 
found in Sect.|6lthat for a constant SFR < 10 yr^^ the 
impact at 30/im was confined to log vL^ < 44.5, which trans- 
lates to X-ray luminosities of log Ljf < 43.5 for the l:l/10pc 
models or log Lx ^ 43.7 for the nozevol/lOpc and 4:l/10pc 
grids. The affected luminosities are low enough that the rest- 
frame 30/im LFs may still be used to discriminate possible 
different evolutions of the AGN type 2/type 1 ratio (Fig. ll5t . 
The method may prove more problematic with galaxies with 
average SFRs > 10 Mq yr^^ 

The sources that produce the majority of the CXRB have 
been identifie d to be p redominantly obscured AGN with 43 < 
log Lx < 44 jBarger et al. 20 05). These objects are typically 
found at z < 1, which may signify a relationship between 
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the obscuring material around an AGN and the star forma- 
tion in the host galax y (Fabian et al. 1998; Franceschini et al. 
ll999tlBallantvne et alji20 06). It is this possibility that makes 
measuring the evolution of R an important goal in the study of 
galaxy evolution. However, there may be a potential difficulty 
in employing the mid-lR LFs calculated here because, as al- 
ready noted, the important luminosity range to distinguish the 
various evolutions is logLx ^ 44, larger than that of the 
objects that dominate the X-ray background. The reason for 
this is that the three evolutions considered have very differ- 
ent predictions for the ratio of type 2/type 1 quasars (Fig.Q, 
but quasars probably have a very different formation and evo- 
lution history than the lower luminosity Seyferts that com- 
prise the bulk of the CXRB. Quasars are more commonly 
found at z > 1, are likely formed by the mergers of mas- 
sive galaxies, and simulations show their host galaxies have 
little ongoing star fo rmation within the galaxy at late times 
(e.g. .'Hopkins et aLl2005) . The evolutions of R described by 
eqns.|3land|4]assume that the AGN covering factor increases 
with z independent of the AGN luminosity. If the attenuat- 
ing gas and dust is connected to the host galaxy SFR then 
this may not be correct at large Lx, and indeed very tenta- 
tive evidence for a slower increase of R with z at larger Lx 
was found by Ballantvne et al. (2006). If this is the case, then 
evidence for this type of luminosity-dependent redshift evo- 
lution will appear in the 30/im LFs, and can be measured us- 
ing the techniques presented here. Furthermore, the effects of 
our luminosity-independent assumption are most important at 
high 2, but Fig.^^show that significant variations in the rest- 
frame 30/im LF due to the different evolutions of R still occur 
at low redshifts. Thus, a determination of the LF at redshifts 
z ^ 0-0.2 will result in an important constraint on the lo- 
cal value of R. Then, tracking the 30/im LF with redshift 
will yield a measurement of how the AGN type 2/type 1 ratio 
changes with z and Lx, a key parameter in understanding the 
formation and evolution of these galaxies. 

7.4. Testing the Unified Model 

To this point, all the calculations presented in this paper 
assumed that the unified model of AGN was correct. That is, 
the fraction of type 2 AGN at a given Lx and z is equal to the 
covering factor of the absorbing material around the AGN. 
These models can also be used to test whether the unification 
scenario is valid at all luminosities. 

As an example, we searched for IRAS objects satisfy- 
ing |6| > 35 degrees, /iy(60/im)> 0.2 Jy and with red- 
shifts between 0.8 and 1.2. This results in 13 objects in 
which we identified 8 quasars using the literature. Two of 
them ([HB89]0235h-164 and [HB89]1308h-326) are blazars 
and are omitted from further consideration. An additional 3 
(IRAS F09121H-2430, PG1248H-401 and B31340H-407) were 
removed because they have weak detections at 60/im and 
the association with the quasars therefore is not firmly es- 
tablished. The remaining three quasars (2MASSiJ01 19560- 
201022, 2MASSiJ1543519H-162422 and PG1206H-459) have 
strong detections by IRAS and the optical data show them to 
have normal QSO SEDs. Assuming z ^ 1, the rest-frame 
30/im luminosities of these quasars are vLi, ^ 46.8. Fig. 1141 
shows that the l:l/10pc evolutionary model, as a result of its 
~3:1 type 2/type 1 ratio (Fig. 0, can reach a 30/im lumi- 
nosity of vLi, = 46.5 at z = 1, but this occurs only when 
the input X-ray luminosity is \ogLx — 48. Yet, the 2- 
10 keV X-ray luminos ity of PG1206H-459 is ~ IQ-^^ erg s^^ 
dPiconcelli et al.l2005h . Indeed, using the three IRAS sources 



to estimate the 30/im LF at vL,^ ~ 46.8 results in a lower- 
limit of > 10^ ^'^ Mpc^"^, about 2 orders of magnitude larger 
than the predicted LF from Fig.[2] The best way of reducing 
this discrepancy would be if the models could produce these 
high IR luminosities at lower AGN input luminosities. 

It seems that the assumed geometric configuration of a rel- 
atively compact absorbing medium situated at 1 pc or 10 pc 
from the AGN cannot produce the necessary 30/im emission 
at high Lx to account for the IRAS limit. Recall that because 
of the assumed density of 10^ cm^'^, the pathlength through 
the absorber is only 15 pc for log A'h = 23 and log Lx = 48 
(taking into account the ionization effects), therefore at high 
AGN luminosities there is only a small column of dust cool 
enough to emit at 30/im. Material at greater distances is 
needed to increase the 30/im emission. A simple way to test 
this idea is to calculate a new grid of Cloudy models with the 
inner radius remaining at 10 pc but the density now lowered 
to 10^ cm^'^. This greatly increases the pathlength of the ab- 
sorbin g m aterial (it is now over 1 kpc in the example above). 
FigureTTTlplots the rest-frame 30/im LFs computed using the 
models with a lower density absorber as well as the previous 
derived ones that assume the more compact geometry. Only 
results for log lyL^ > 44.7 are plotted for the low density ab- 
sorber. The greater pathlength afforded by the lower density 
cloud now produces enough 30/im emission to account for 
the IRAS lower-limit. The right-hand panel of Fig.ll7lshows 
the same results but for the rest- frame 8/im LF. At this wave- 
length, the differences in the LF between the two assumed 
densities are small. 

The implication of this simple test is that the high luminos- 
ity quasars require absorption and remission of nuclear radia- 
tion over a large range of radii, while the current data available 
for the lower luminosity quasars and Seyfert galaxies can be 
adequately described by a more compact absorber. The type 
2/type 1 ratio at the IRAS luminosities for the three evolu- 
tions considered in this paper varies from ^1:1-4:1 (Fig.Q, 
a range that is consistent with the latest observational con- 
straints dMartmez-Sansigre et aL.2005-.2006) . Under the uni- 
fication model, this translates into covering fractions of 0.5- 
0.8. Producing high covering fractions is significantly simpler 
if the obscuring material is close to the AGN, and not spread 
out over a large range of radii, as implied by the above results. 
Clearly, the change in mid-IR LF with luminosity (in addition 
to the possible redshift evolution concentrated on in previous 
sections) contains fundamental clues on the nature and geom- 
etry of the obscuring medium and therefore the host galaxy. 
Future observations are needed for this technique to fulfill this 
promise. 

We conclude that the unification model based on a relatively 
compact absorber is consistent with the available data for 
moderate AGN luminosities (Seyfert galaxies) but fails when 
extending this to the high luminosity quasars. This adds to 
the evidence that the evolution and fueling of quasars are very 
different from those that drive the lower luminosity Seyfert 
galaxies. 

8. SUMMARY 

The principle conclusion of this work is that a measure- 
ment of the redshift evolution of the mid-IR AGN LF will 
determine to what extent the AGN type 2/type 1 ratio, R, 
is evolving with z. The specific predictions shown above 
were of most use for the existing Spitzer MIPS bands, but are 
equally valid for future observations with JWST. Predictions 
were shown for three different evolutions of R, but the current 
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data were unable to distinguish among them. 

The calculations performed here were relatively novel, as 
they employed the photoionization code Cloudy to directly 
connect the X-ray properties of AGN to the IR. This tech- 
nique, although not appropriate for fitting SEDs of individual 
objects, yielded A^n-averaged SEDs that well described the 
average properties of AGN found in the deep surveys. In addi- 
tion, varying the parameters of the models allows constraints 
to be placed on the physical conditions for an average AGN 
absorber. For example, comparing results assuming two dif- 
ferent distances for the average inner radius of the attenuating 
gas and dust (1 pc and 10 pc), yielded compelling evidence 
that the rest frame mid-lR properties of a large ensemble of 
AGN were best described when the inner radius of the ob- 
scuring material was ^ 10 pc from the AGN. Future work 
will cover a greater parameter space and split the data into 
both redshift and luminosity bins. In this way, properties of 
average AGN can be derived and followed with minimal in- 
terference from the intrinsic object-to-object variability. 

Finally, it is worth emphasizing the importance of determin- 
ing the evolution of the AGN type 2/type 1 ratio, as it may 
be directly related to processes within the host galaxy. Be- 
tween z ^ and 1, most AGN are relatively low-luminosity 
Seyfert galaxies and therefore such studies probe a differ- 



ent mode of galaxy formation and evolution than the higher- 
luminosity quasars at high redshifts. Evidence for such a 
difference was uncovered by estimating the rest-frame 30/im 
LF at high quasar luminosities. This LF was more than two 
orders of magnitude above the predictions, which assumed 
that a compact unification model was valid at all luminosities. 
To account for this observational limit absorbing material at 
much larger radii was invoked, suggesting a different mode of 
obscuration. Seyfert galaxies, in contrast, can still be success- 
fully described by the compact unified model. 
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Fig. 1 . — (a) Contours of the AGN type 2/type 1 ratio (R) as a function of Lx and redsliift for the evolution described by eqn.|3] Adapted from'Ballantvne et al| 
120061) . Recall that the z-evolution is halted at ^ = 1. (b) Same as (a), but shows contours of R for the more rapid evolution of eqn.|4] (c) As in (a) and (b), but 
for a model with no redshift evolution (Eg.lsl. 
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Fig. 2. — (Top) Example of a A^H-averaged SED. In this case, the illuminating AGN has a 2-10 keV luminosity of 10^^ erg s~^, and the inner radius of the 
absorbing gas is 1 pc from the ionizing source. This spectrum is taken from the 4:1 evolutionary grid (eqn.|3) when z = 0.45. Thus the fraction of type 2 AGN 
is /2 = 0.754. The solid black line is the final A'jj -averaged SED when the simple Nn distribution is assumed, while the dotted black line is the result when 
the Risaliti et al. Nn distribution is used. The Risaliti et al. distribution increases the emission longward of ~30^m, and diminishes it between ~ 1 and 30/im, 
as compared to the result with the uniform A'^h distribution. (Bottom) As above, but with r = 10 pc. 
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Fig. 3. — (a) to 2-8 keV flux ratio versus X-ray hardness for afl 4:l/lpc models with ^ < 3 and logLx < 48. The fluxes are taken from the 

unifled SEDs, i.e., before the flnal average over the A^h distribution. The blue region covers models with log A'^h > 24, the red area shows the results when 
23 < log A'^H < 24, and the black region is derived from models with log N^i < 23. The green squares are data from the Chandra Deep Field South taken 
from Fig. 2 in the paper by Rigbv et al. 1 2004). (b) As in (a) except only models with log Lx < 43 are plotted, (c) As in (a) except the results are taken from 
the 4:l/10pc grid, (d) As in (b) except the results are taken from the 4:l/10pc series of models. 
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Fig. 4. — As in Figurelsja) except the results for the final A^n-averaged SEDs are covered by the black region. As a result of the averaging, the final SEDs no 
longer span a large range of X-ray hardness. The red area coiresponds to AGN with z < 1 and log Lx < 44. These sources produce the largest contribution to 
the CXRB and so must have a hardness ratio corresponding to a F = 1.4 power-law (i.e., log(/2_8 kcV / /o.5-2 kcv) = 0.36, plotted as a dotted line). Sources 
at higher redshift will appear softer. 
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log 

Fig. 5. — The ratio of input Lx to vL^(5.1 fim) plotted versus log A'^h using the unified SEDs from the 4:l/lpc (top) and 4:l/10pc (bottom) grids. The 
points are color-coded by Lx- black triangles denote 46 < logLx < 48, red triangles denote 44 < logLx < 46, and blue triangles show results for 
41.5 < log Lx < 44. Only models at z = 0.7 are plotted, but as rest-frame luminosities are being plotted, there is little change at other redshifts. The green 
squares are data from both type 1 and 2 AGN taken from Rigbv et al, ( 2006), while the cyan stars are data from Lutz et al. 1 20()4). All of the observed data use 
6/xm luminosities. 
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Fig. 6. — A color-color diagram for the 4 IRAC bands from Spitzer, after lLacv et'al] l2004h . (Left) The black crosses are the colors predicted from the 4:l/lpc 
iVn -avera ged SEDs when ^ < 3. All Lx are included in this diagram. The region enclosed by the three dashed lines is the AGN selection region used by 
ILacv et alj t20Q41 (Right) Same as left-hand panel, but results are plotted for the 4:l/10pc grid. 




Fig. 7. — A color-color diagram across the MIPS-IRAC bands from Spitzer, afterl Lacv etaH 1200 411. (Left) The black crosses are the colors predicted from the 
L1/1£c_7Vh -averaged SEDs when z < 3. All values of Lx are included in the figure. The red squares are the observed colors of a variety of AGN plotted by 
ILacv et alj 120041 in his Figure 2. (Right) Same as left-hand panel, but results are plotted for l:l/10pc grid. 
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Fig. 8. — Histogram of the 24/im/8/im flux ratio for sources in the nozevol/lpc (black solid line), the nozevol/lOpc (black dashed line), the l:l/lpc (red 
solid line ) and the l:l/10pc (red dashed line) grids. The dotted-Hne plots the histogram of flux ratios observed from X-ray selected AGN in the XBootes survey 
IBrand et aIji20Q6i) . Only the models with /0.5— 7 kcV > 7.8 X 10~^^ erg cm~^ s~^, the flux limit of the XBootes survey, and /^(24/im) > 0.3 mjy are 
plotted in the figure. 
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Fig. 9. — (Top) The predicted 8^m flux is plotted against the 0.5-10 keV X-ray flux. The black crosses denote results from the 4:l/lpc grid, while the blue 
triangles are from the 4:l/10pc grid. Predictions were only plotted for z < 3. The red squares are the observed fluxes from X-ray sources in the Extended Groth 
Strip taken from Fig. 7 in the paper bv lBarmbv et alji2006D . (Bottom) As above, except the 3.6/im flux is plotted against the 0.5-10 keV flux. 
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Fig. 10. — The predicted 2-10 keV flux is plotted against the 24/im flux. The symbols are the same as Fig.|9] except the red squares are data from the Chandra 
Deep Field South taken from Fig. 1 in the paper by Alonso-Henero et al. 1 2006 ). The hatched regions denote the extrapolation of the hard X-ray to mid-IR ratios 
for low-redshift X-ray selected AGN ('Piccionotti AGN') and local starburst galaxies (see iAlonso-Herrero et alj2004.i2006l) . 
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Fig. 1 1 . — Cumulative number count distributions for 5.7, 8, 24 and IQfim. The black lines plot the predicted distributions from the models with the inner 
radius of the attenuating gas and dust at r = 1 pc from the ionizing source, while the red lines show the results when r = 10 pc. The solid lines denote the no z 
evolution model, while the dotted and dashed lines plot the 1 : 1 and 4: 1 evolving models respectively. Green dashed lines show results for the non-evolving model 
computed with ISM grains. They are only relevant for the 5.7 and 8/im panels. The simple A'^jj distribution has been used in constructing the distributions. The 
data points in the 5.7, 8, and 24/xm panels are taken from the GOODS survey of X-ray selected AGN iiTreister et aU2006ft . 
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Fig. 12. — Input hard X-ray luminosity versus predicted IR luminosity for the Afjj-averaged SEDs in the nozevol/lpc (left panel) and nozevol/lOpc (right 
panel) grids. As indicated in the figures, the various line styles denote different IR wavelengths. The range of luminosity spanned in the infrared is smaller than 
the range in the X-ray, with the longer wavelengths covering smaller ranges. 




Fig. 13. — (Left) Predicted rest-frame AGN 8/xm luminosity functions for the three different evolutions of the type 2/type 1 ratio. The black lines show 
the results for the nozevol/lpc grid, the red lines denote the l:l/lpc grid, and the blue lines plot the 4:l/lpc results. Three different redshifts are plotted using 
different line styles. The solid, dotted and dashed lines show the luminosity functions at z = 0.2, 0.6 and 2.0, respectively. The lower panel plots the ratio of 
the evolving luminosity functions to the non-evolving one for the three different redshifts. The colors denote the different evolving models as in the upper panel. 
For example, in the lower-panel the red-dashed line plots the ratio of the predicted l:l/lpc luminosity function at 2 = 2 t o the p redicted nozevol/lpc luminosity 
function at z = 2. The data points are the measured luminosity function of type 1 quasars determined by Brown et^H 12006^ . with the stars denoting objects 
with 1 < 2 < 5 and the diamonds are objects with 1.5 < z < 2.5. (Right) As in the left-hand side, but for models with the inner radius of the absorbing gas 
and dust at a distance of 10 pc from the AGN. 
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Fig. 15. — The predicted rest-frame 30um luminosity functions from the 1:1 grids at z = 0.6. The solid line has no contribution from star-formation and 
is identical to the equivalent curve in Fig. 1141 The other line styles show the prediction when a constant SFR is included using the dust emission spectrum 
of eqn.Q dotted-line=0.5 Mq yr~^, short-dashed line=l Mq yr~^, long-dashed line=5 Mq yr~^, short dot-dashed line=10 Mq yr^^, long dot-dashed 
line=100 Mq yr~^ . The left(right)-hand panel contains the results when the inner radius of the attenuating material is at 1(10) pc from the AGN. 
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Fig. 17. — (Left) Plot of the rest-frame 3()/im LF at z = 0.2, 0.6 and 1.0 (line styles and low-luminosity data point as in 

Fig.ini- The LFs were computed 

using the nozevol/lOpc model grid. The arrow denotes the lower-l imit to the high-luminosity rest-frame 30^m LF at z = 1 derived from IRAS sources (see t ext). 
The black lines are identical to the coiresponding ones from Fig. ll4l and plot the predictions using the assumed compact geometry described in Sect. 13.21 The 
red lines show the high-vL^ results when the density of the absorbing cloud is lowered to 100 cm^'^. The absorber extends to much greater distances than the 
previous case, and the refo re enhances the emission in the mid-IR. (Right) As in the the left-hand plot, but for the rest-frame 8/im LF. The line styles and points 
are the same as in Fig. 1131 



